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ABSTRACT 

We study the "minimal" cooling scenario of superfluid neutron stars with 
nucleon cores, where the direct Urea process is forbidden and the enhanced 
cooling is produced by the neutrino emission due to Cooper pairing of neu- 
trons. Extending our previous consideration (Gusakov et al. 2004a), we include 
the effects of accreted envelopes of light elements. We employ phenomenologi- 
cal density-dependent critical temperatures T cp (p) and T cnt (p) of singlet-state 
proton and triplet-state neutron pairing in a stellar core, as well as the criti- 
cal temperature T cns (p) of singlet-state neutron pairing in a stellar crust. We 
show that the presence of accreted envelopes simplifies the interpretation of 
observations of thermal radiation from isolated neutron stars in the scenario of 
Gusakov et al. (2004a) and widens the class of models for nucleon superfluidity 
in neutron star interiors consistent with the observations. 
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1 INTRODUCTION 

New observations of thermal radiation from isolated middle-aged neutron stars (e.g., Pavlov, 
Zavlin & Sanwal 2002, Pavlov & Zavlin 2003) initiated further development of the cooling 
theory of these objects. Its main aim is to interpret the data and constrain still poorly known 
properties of dense matter in neutron star cores, such as the composition, the equation 

* E-mail: kam@astro.ioffe.ru (ADK), gusakov@astro.ioffe.ru (MEG), yak@astro.ioffe.ru (DGY), ognedin@astronomy.ohio- 
state.edu (OYG) 



2 A. D. Kaminker, M. E. Gusakov, D. G. Yakovlev, O. Y. Gnedin 

of state and nucleoli superfluidity (e.g., Yakovlev & Pethick 2004, Page et al. 2004 and 
references therein). 

It is well-known (e.g., Yakovlev & Pethick 2004) that theoretical models of non-superfluid 
neutron stars which possess nucleon cores and cool via the modified Urea process of neutrino 
emission cannot explain the observations. Some neutron stars (e.g., RX J0822-4300 and 
PSR B 1055-52) are much warmer than predicted by these theories, while others (e.g., the 
Vela pulsar or the compact source in CTA 1) are much colder. Warmest objects can be 
treated as relatively low-mass neutron stars with strong proton (e.g., Kaminker, Haensel & 
Yakovlev 2001) or neutron (e.g., Gusakov et al. 2004b) pairing in their cores. Strong pairing 
suppresses the modified Urea process and makes the stars warmer. Coldest stars should have 
higher neutrino emission than the emission provided by the modified Urea process. They 
are usually treated as massive neutron stars which cool either via the powerful direct Urea 
process in nucleon (or nucleon/hyperon) matter or via similar processes in kaon-condensed, 
pion-condensed, or quark matter in their inner cores. 

Recently Page et al. (2004) and Gusakov et al. (2004a) proposed new scenarios of neutron 
star cooling which involve only standard physics of neutron star interiors. The neutron star 
cores are assumed to contain nucleons (no exotic forms of matter) with the forbidden direct 
Urea process. Some enhancement of the cooling can be provided by neutrino emission due 
to Cooper pairing of nucleons. Page et al. (2004) called their cooling scenario the "Minimal 
Cooling Model" (for its simplicity). We will also use this very properly chosen name for the 
scenario of Gusakov et al. (2004a) that is based on the same assumptions (but differs in 
their realization; see below). 

According to our previous paper (Gusakov et al. 2004a) the enhanced cooling is produced 
by the neutrino emission due to Cooper pairing of neutrons in the cores of massive neutron 
stars, while warmest objects are thought to be low-mass stars with strong proton pairing in 
their cores. We assumed a phenomenological model of strong density-dependent singlet-state 
proton pairing with the critical temperature T cp (p) that has the maximum value T™ ax > 5.0 x 
10 9 K. We also assumed a phenomenological model of moderate triplet-state neutron pairing 
T cnt (p) with the maximum critical temperature T™ t ax ~ 6.0 x 10 s K shifted to higher p, where 
proton pairing dies out. We were able to interpret all the data but under stringent constraints 
on the density dependence of T cnt (p). 

The present paper extends our previous analysis. We use the same equation of state of 
matter in neutron star interiors (Douchin & Haensel 2001) and the same model of triplet- 
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state neutron pairing. However, in addition, we take into account the effects of surface 
layers of light (accreted) elements (H and/or He), as well as singlet-state neutron pairing 
Tcns(p) in the stellar crust. The effects of accreted envelopes allow us to lower proton pairing 
(Tcp iax ~ 10 9 K) required to explain the data. This weaker proton pairing is consistent 
with recent microscopic calculations of proton critical temperatures by Zuo et al. (2004) 
and Takatsuka & Tamagaki (2004) (although some other calculations predict much stronger 
proton pairing; e.g., Lombardo & Schulze 2001; also see references in Yakovlev, Levenfish & 
Shibanov 1999, and a recent paper by Tanigawa, Matsuzaki & Chiba 2004). 

Let us emphasize the difference of cooling scenarios of Page et al. (2004) and Gusakov 
et al. (2004a). In particular, Page et al. (2004) used several selected models of triplet- 
state neutron pairing provided by microscopic theories. Corresponding cooling curves do 
not depend sensitively on neutron star mass and do not allow the authors to explain all the 
data in the frame of one physical model of neutron star interiors. In contrast, Gusakov et al. 
(2004a) used phenomenological models of triplet-state pairing and succeeded to explain all 
the data (although under stringent constraints on these models; see their paper for details). 
Note that Page el al. (2004) analyzed the effect of accreted envelopes on their minimal 
cooling models but our models are different and require separate analysis. Our main aim 
is to interpret all the data assuming the same physics (equation of state and superfluid 
properties) in the interiors of all neutron stars. 



2 OBSERVATIONS 

Table 1 summarizes observations of isolated neutron stars, whose thermal surface radia- 
tion has been detected or constrained. We present the estimated stellar age t, the effective 
surface temperature T s °° and the surface thermal luminosities (as detected by a dis- 
tant observer). The data on t and T s °° are described by Gusakov et al. (2004a) in more 
detail, with two exceptions. First, following Slane et al. (2004a), we slightly lower the upper 
limit on the surface temperature T s °° of PSR J0205+6449 in the supernova remnant 3C 58 
(T s °° < 1.02 MK instead of 1.1 MK). Second, we include into consideration the central X-ray 
source RX J0007. 0+7303 in the supernova remnant CTA 1. 

For PSR J0205+6449 we adopt the age of the historical supernova SN 1181 (t &820 
yr). However notice, that recently Chevalier (2004, 2005) presented arguments in favor for 
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Table 1. Observational limits on surface temperatures and thermal luminosities of isolated neutron stars 



Source 


t [kyr] 


T s °° [MK] 


Confid. 


Refs. c ) 


l gL oo [crg/s] 


PSR J0205+6449 (in 3C 58) 


0.82 


<1.02 b ) 


99.8% 


S04a 


< 33.29 


PSR B0531+21 (Crab) 


1 


<2.0 b) 


99.8% 


W04 


< 34.45 


RX J0822-4300 


2-5 


1.6-1.9 a > 


90% 


ZTP99 


33.9-34.2 


IE 1207.4-5209 


3-20 


1.4-1.9 a ) 


90% 


ZPS04 


33.67-34.20 


RX J0007.0+7303 (in CTA 1) 


10-30 


< 0.66 ft > 




H04 


< 32.54 


PSR B0833-45 (Vela) 


11-25 


0.65-0.71 a > 


68% 


P01 


32.19-32.67 


PSR B 1706-44 


~17 


u -° -0.34 


68% 


M04 


31.66-32.94 


PSR J0538+2817 


30 ±4 


~ 0.87 a > 




ZP04 


32.32-33.33 


PSR B0633+1748 (Geminga) 


~340 


~ 0.5 b > 




K05 


31.34-32.37 


RX J1856.4-3754 


~500 


<0.65 




G04a 


< 32.5 


PSR B1055-52 


~540 


~0.75 b ) 




PZ03 


32.05-33.08 


RX J0720.4-3125 


~ 1300 


~ 0.51 




MZH03 


31.37-32.40 



a ) Inferred using a hydrogen atmosphere model 

Inferred using the black-body spectrum 
c > S04a - Slane et al. (2004a); W04 - Weisskopf et al. (2004); ZTP99 - Zavlin, Triimper & Pavlov (1999); ZPS04 - Zavlin, 
Pavlov & Sanwal (2004); H04 - Halpern et al. (2004); P01 - Pavlov ct al. (2001); M04 - McGowan et al. (2004); ZP04 - 
Zavlin & Pavlov (2004); K05 - Kargaltsev ct al. (2005); G04a - see Gusakov ct al. (2004a); PZ03 - Pavlov & Zavlin (2003); 
MZH03 - Motch, Zavlin & Haberl (2003) 

a larger age of the pulsar wind nebula in 3C 58 (t =2400±500 yr). Were this the actual age 
of the neutron star, its interpretation would be easier. 

For the source RX J0007. 0+7303 we adopt the age of its host supernova remnant CTA 1 
(Gl 19.5+10.2). According to Slane et al. (2004b), the age is t = 13 kyr. Following Halpern et 
al. (2004) we assume the neutron star age limits 10 kyr < t < 30 kyr. As for RX J0205+6449, 
the Crab pulsar and RX J0007. 0+7303, no thermal radiation component has been detected 
from these objects, and only the upper limits on T s °° have been set (Slane et al. 2004a, 
Weisskopf et al. 2004, Slane et al. 2004b, Halpern et al. 2004). 

The surface temperatures of some sources from Table 1 (labeled by ")) have been obtained 
by fitting their thermal radiation spectra with hydrogen atmosphere models. Such models are 
more consistent with other information on these sources (e.g., Pavlov & Zavlin 2003) than 
the blackbody model. For other sources (e.g., for the Geminga pulsar and PSR B1055-52, 
labeled by we present the values of T s °° inferred using the blackbody spectrum because this 
spectrum is more consistent for these sources. The surface temperature of RX J1856. 4-3754 
is still uncertain. Following Gusakov et al. (2004a) we adopt the upper limit T s °° < 0.65 MK. 
Finally, T s °° for RX J0720.4-3125 is taken from Motch et al. (2003), who interpreted the 
observed spectrum with a model of a hydrogen atmosphere of finite depth. Note also the 
new results by Kargaltsev et al. (2005) for Geminga presented in Table 1. These authors 
confirm the observational value of T s °° reported by Zavlin & Pavlov (2004). Taking into 
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account systematic uncertainties of T s °° discussed by Kargaltsev et al. (2005) we retain 20% 
errorbars adopted by Gusakov et al. (2004a) and erroneously referred to 90% confidence level 
in their Table 1. Following Gusakov et al. (2004a), the same 20% errorbars will be adopted 
for PSR J0538+2817, PSR B1055-52, and RX J0720.4-3128. 

As noted by several authors (e.g., Page et al. 2004), it may be instructive to compare 
the cooling theory with measured values of stellar thermal surface luminosities L™, rather 
than with T s °°. The data on are also collected in Table 1. The luminosity is related to 
the effective surface temperature via 

L? =4no-Rl(Tn\ (1) 

where a is the Stefan-Boltzmann constant, Roo = Rj \J\ — 2GM/(c 2 R) is the so called 
apparent radius of a neutron star (as would be detected by a distant observer if a telescope 
could resolve the star), R is the circumferential radius, and M the gravitational stellar mass. 
Thus, the luminosity is determined by the effective temperature and neutron star radius; an 
uncertainty in is produced by uncertainties in T s °° and R^. We have already described 
the values of T s °°. As for the values of R^, we vary them (with two exceptions indicated 
below) within the reasonable theoretical interval for neutron star radii, -Roo=ll-16 km; while 
translating R into R^ we always set M = 1.4 M Q . 

In Table 1 the upper limits on for PSR J0205+6449, the Crab pulsar, RX J0007.0+7303, 
and RX J1856. 4-3754 were obtained assuming R^ = 16 km. 

The luminosities of RX J0822-4300 and IE 1207.4-5209 have been calculated from the 
values of T s °° obtained by Zavlin et al. (1999) and Zavlin et al. (2004), respectively. We have 
taken the same fixed radius R = 10 km (R^ = 13 km) which was used by the cited authors 
to fit the observed spectra with the hydrogen atmosphere models. All other values of in 
Table 1 have been obtained by varying R^ within the interval R^ =11-16 km. 

The central values of have been calculated taking into account the central values of 
T s °° from Table 1 and the values of R (or R^) obtained in cited papers from spectral fits, 
except for the Vela pulsar, where we set R^ = 13 km. For PSR B1706-44, PSR J0538+2817, 
and RX J0720. 4-3125 these values of R have been taking 12 km, 10.5 km, and 10 km, as 
suggested by McGowan et al. (2004), Zavlin & Pavlov (2004), and Motch et al. (2003), 
respectively. For the Geminga pulsar we have used the value R = 10.6 km from Zavlin & 
Pavlov (2004), and for PSR B 1055-52 we set R = 13 km from Pavlov & Zavlin (2003). 

In all the cases, the limits of presented in Table 1 seem to be rather uncertain. 



6 A. D. Kaminker, M. E. Gusakov, D. G. Yakovlev, O. Y. Gnedin 

Although, in principle, the luminosities can be measured/constrained more accurately 
than T s °° (by exact measuring the distance and the bolometric thermal flux), it is not so for 
the sources collected in Table 1 mainly due to large uncertainties in measured distances to 
the sources (see, e.g., Page et al. 2004). Nevertheless, comparing observed and theoretical 
luminosities of cooling neutron stars seems to be useful. Our limits of are in reasonable 
agreement with corresponding limits given by Page et al. (2004). The main differences refer 
to the Geminga pulsar and IE 1207.4-5209. In the first case the limits of Lf presented 
by Page et al. (2004) correlate with too low apparent radius of the star, i?^ < 6 km, 
for the temperature limits adopted in their paper. In the second case Page et al. used 
L s °° = L s {R/Roo) 2 with L s = 5.0t 4 j x 10 33 erg s~ l from Zavlin, Pavlov & Trumper (1998). 
The value of L s was possibly underestimated by Zavlin et al. (1998), because their value of 
T s was indicated later by Zavlin et al. (2004) as T s °°. 

Also, let us note that the radii of our neutron star models used for the cooling calculations 
presented below are consistent with the radii used for the interpretation of the data. 

3 PHYSICS INPUT AND CALCULATIONS 

The cooling calculations have been done using our general relativistic cooling code (Gnedin, 
Yakovlev & Potekhin 2001). At the initial cooling stage (t < 100 yr) the main cooling 
mechanism is the neutrino emission but the stellar interior stays highly non- isothermal. At 
the next stage (10 2 yr < t < 10 5 yr) the neutrino emission is dominant but the stellar 
interior is isothermal. Later (t > 10 5 yr) the star cools predominantly through the surface 
photon emission. 

Following Gusakov el al. (2004a) we adopt the moderately stiff equation of state for the 
neutron star matter suggested by Douchin & Haensel (2001). In this case a neutron star 
core (a region of density p > 1.3 x 10 14 g cm -3 ) consists of neutrons with the admixture 
of protons, electrons and muons. All constituents exist everywhere in the core, except for 
muons which appear at p > 2.03 x 10 14 g cm~ 3 . The most massive stable star has the 
(gravitational) mass M = M max = 2.05 M , the central density p c = 2.9 x 10 15 g cm -3 , and 
the (circumferential) radius R = 9.99 km. The parameters of neutron stars with some other 
masses are given by Gusakov et al. (2004a). 

The employed equation of state forbids the powerful direct Urea process of neutrino 
emission (Lattimer et al. 1991) in all stable neutron stars (M ^ M max ). Accordingly, a non- 
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superfluid neutron star of any mass in the range M & < M ^ M max (without any accreted 
envelope) will have almost the same (universal) cooling curve T s °°(i) (the dotted curve in the 
right panel of Fig.QJ. At the neutrino cooling stage, this curve is determined by the modified 
Urea process and is almost independent of the equation of state in the stellar core (see, e.g., 
Yakovlev & Pethick 2004 and references therein). As seen from Fig. Q this universal cooling 
curve cannot explain the data. We will show that all the data can be explained assuming 
nucleon superfluidity in the internal layers of neutron stars and the presence of accreted 
envelopes (of light elements). 

Following the standard procedure (Gudmundsson, Pethick & Epstein 1983) our code 
calculates heat transport in the neutron-star interior (p > pb — 10 10 g cm" 3 ) and uses the 
predetermined relation between the effective surface temperature T s and the temperature 
7b at the bottom of the surface heat-blanketing envelope (p < pb). We use the relation 
calculated by Potekhin, Chabrier, & Yakovlev (1997) and updated by Potekhin et al. (2003). 
We will employ the models of blanketing envelopes made of iron (which is the standard 
assumption) and envelopes containing light elements. 

The detailed description of these models is given by Potekhin et al. (2003). The thermal 
energy in the heat-blanketing envelope is mainly conducted by electrons. The thermal con- 
ductivity of electrons which scatter off lighter ions in the accreted envelope is higher than 
the conductivity in the iron envelope. This means that the accreted envelope is more heat 
transparent than the iron one, resulting in higher T s for the same Tb. This rise of the surface 
temperature depends on Tb and AM, the mass of light elements (hydrogen and/or helium, 
with a possible carbon/oxygen layer at the bottom of the accreted envelope as a result of 
nuclear burning of lighter elements). Potekhin et al. (1997, 2003) varied the boundaries of 
layers containing different elements within physically reasonable limits and found that the 
resulting relation between T s and Tb is remarkably insensitive to these variations and de- 
pends mainly on AM. However, AM cannot exceed ~ 10~ 7 M, because at higher AM the 
bottom density of the accreted envelope would exceed 10 10 g cm~ 3 . At such high densities, 
light elements (including carbon/oxygen) would rapidly transform into heavier ones. 

At the neutrino cooling stage Tb is governed by the neutrino emission from the stellar 
interior and is almost independent of conductive properties in the heat-blanketing envelope. 
In contrast, at the photon cooling stage the star with the accreted envelope has lower Tb 
and, consequently, lower T s due to higher heat transparency of the surface layers. This leads 
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to faster photon cooling through the surface (for not too cold stars; see, e.g., Potekhin et al. 
1997). 

The cooling of a neutron star is sensitive to superfluidity of nucleons in the stellar core and 
to superfluidity of free neutrons in the inner stellar crust. Any superfluidity is characterized 
by its own density-dependent critical temperature T c (p). Microscopic theories predict mainly 
(i) singlet-state ( 1 So) pairing of neutrons (T c = T cns ) in the inner crust and the outermost 
core; (ii) 1 So proton pairing in the core (T c = T cp ); and (iii) triplet-state ( 3 P2) neutron 
pairing in the core (T c = T cnt ). These theories give a large scatter of critical temperatures, 
from ~ 10 10 K to ~ 10 8 K and lower, depending on a nucleon-nucleon interaction model and 
a many-body theory employed (e.g., Lombardo & Schulze 2001, Yakovlev et al. 1999; also 
see recent papers by Schwenk & Friman 2004, Takatsuka & Tamagaki 2004, Zuo et al. 2004, 
Tanigawa et al. 2004). Because of these huge theoretical uncertainties, we will not rely on any 
specific microscopic results but will treat T cp (p) and T cn (p) as phenomenological functions 
of p (which can be varied in physically reasonable limits). Our aim will be to constrain these 
functions by comparing theoretical cooling curves with the observations. 

Superfluidity of nucleons affects the heat capacity and suppresses neutrino processes such 
as Urea and nucleon-nucleon bremsstrahlung processes (as reviewed, e.g., by Yakovlev et al. 
1999). It also introduces an additional neutrino emission mechanism associated with Cooper 
pairing of nucleons (Flowers, Ruderman & Sutherland 1976). All these effects of superfluid- 
ity are incorporated into our cooling code. While calculating the neutrino emission due to 
Cooper pairing of protons we use phenomenological values of weak-interaction parameters 
renormalized by many-body effects (the same as in Gusakov et al. 2004b). 

In the left panel of Fig. ^we plot models for nucleon pairing adopted in our calculations: 
one model nsl of strong singlet-state pairing of neutrons (with the peak of T cns (p) approxi- 
mately equal to T™g X pa 7 x 10 9 K); three models of proton pairing - strong pi, moderately 
strong p2, and moderate p3 (T™ ax ~ 6.8 x 10 9 K, 1.5 x 10 9 K, and 7.5 x 10 8 K, respectively); 
and one model ntl of moderate triplet-state neutron pairing (T™ ax ~ 6 x 10 8 K). Models 
pi and ntl are the same as in Gusakov et al. (2004a). Now we add models of weaker proton 
pairing (particularly, p2). Strong proton pairing has been predicted in a number of publi- 
cations (e.g., Tanigawa et al. 2004) while other publications predict much weaker proton 
pairing (e.g., Zuo et al. 2004, Takatsuka & Tamagaki 2004). 

As seen from the right panel of Fig. Q proton pairing p2 combined with strong crustal 
superfluidity of neutrons nsl results in too cold low-mass neutron stars. The neutrino emis- 
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Figure 1. Left: Density dependence of critical temperatures for three models pi, p2, and p3 of singlet-state proton pairing 
and one model ntl of triplet-state neutron pairing in the neutron star core, as well as model nsl of singlet-state neutron pairing 
in the stellar crust. Vertical dot-and-dash lines indicate the crust-core boundary and the central density of the maximum-mass 
neutron star. Right: Observational limits of surface temperatures of neutron stars (Table 1) as compared with theoretical 
cooling curves. Three solid curves refer to neutron stars of different masses (indicated near the curves) with nucleon pairing 
p2, ntl, and nsl. The dot-and-dashed curve refers to an M = I.IUMq neutron star with the same superfluidity, but with the 
accreted envelope of the mass AM = 10 — 8 M. The dotted curve is for a non-superfluid star of the same mass without any 
accreted envelope. 

sion due to Cooper pairing of protons in the core and of neutrons in the inner crust (see 
Section HJ) accelerates the cooling and does not allow us to explain the observations of the 
young and hot neutron stars, RX J0822-4300 and IE 1207.4-5209. However, this cooling 
scenario is consistent with the observations of the old and warm neutron stars, PSR 1055-52 
and RX J0720. 4-3125. Accreted envelopes can rise the surface temperatures of middle-aged 
neutron stars and explain the observations of RX J0822-4300 and IE 1207.4-5209. This is 
demonstrated by the dot-and-dashed cooling curve for the low-mass star with the accreted 
envelope of the mass AM = 10 _8 M. 

Our interpretation of the neutron stars coldest for their age (PSR J0205+6449 in 3C 58, 
RX J0007. 0+7303 in CTA 1, the Vela and Geminga pulsars) remains the same as in Gusakov 
et al. (2004a). These objects can be treated as massive neutron stars (M > 1.9 M©) with 
moderate triplet-state neutron pairing ntl in their inner cores where proton pairing pi (as 
well as p2 and p3) dies out (the left panel of Fig.^). Our phenomeno logical pairing model ntl 
seems specific (shifted to too high densities p). However similar models have been obtained 
from microscopic theories (e.g., see the curve m* = 0.73 in Fig. 1 of Takatsuka & Tamagaki 
1997). 

In this way we come to the same three distinct classes of cooling neutron stars as in 
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Figure 2. Cooling of low-mass (1.111 Mq) neutron stars with proton pairing (models pi, p2 or p3) and neutron pairing ntl 
(weak at low densities, Fig. Q in the stellar core (thin short-dashed, solid, and long-dashed curves, respectively), with the aim 
to interpret the observations of neutron stars hottest for their age (RX J0822-4300, IE 1207.4-5209, PSR B1055-52, and RX 
J0720.4— 3125). Thick curves include, in addition, the effects of crustal neutron superfluidity nsl. Left: No accreted envelopes. 
Right: Thick curves are the same as in the left panel. Thin dot-and-dashed and solid curves include, in addition, the effects of 
accreted envelope of the mass AM = 1Q~ S M. 

Gusakov et al. (2004a) (and generally as in Kaminker et al. 2002). The first class contains low- 
mass stars whose surface layers are composed either of iron or of light elements (solid or dot- 
and-dashed cooling curves, respectively, for the M = 1.111 M Q star in Fig. [J). Another class 
contains high-mass stars which show enhanced cooling (the solid curve for the M = 1.994 M Q 
star) produced by the neutrino emission due to Cooper pairing of neutrons. Finally, there is 
a class of medium-mass neutron stars (the solid curve for the M = 1.37 M Q star) which show 
intermediate cooling. Their cooling curves fill in the space between the upper curve for low- 
mass stars and the lower curve for high-mass stars. These curves explain the observations 
of PSR B1706-44, PSR J0538+2817, and RX J1856.4-3754. 

4 COOLING OF LOW-MASS NEUTRON STARS 

As has been shown in Section 01 the presence of light elements on the surfaces of the younger 
and hotter neutron stars, RX J0822-4300 and IE 1207.4-5209, can allow us to explain their 
observations if we assume moderately strong proton pairing p2 in their interiors. This pairing 
is also consistent with the observations of the old and warmest sources, PSR B1055-52 and 
RX J0720. 4-3125. We interpret all these sources as low-mass neutron stars. Let us analyze 
the main cooling regulators of such stars. 
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In our case, triplet-state neutron pairing in low-mass stars is weak. For the adopted 
equation of state of Douchin & Haensel (2001), this implies T cn (p) < 2 x 10 8 K at p < 8 x 10 14 
g cm -3 . Under this condition, neutron pairing does not affect the cooling of low- mass stars 
(M < 1.1 M©) at least at the neutrino cooling stage. The thin short-dash line in the left 
panel of Fig. |21 shows that (in the absence of crustal pairing) strong proton pairing pi is 
needed to explain the data on all neutron stars hottest for their age (Gusakov et al. 2004a). 
In contrast, cooling curves for moderately strong proton pairing p2 (the thin solid line) 
and moderate pairing p3 (the thin long-dashed line) go essentially lower than the curve for 
pairing pi, being inconsistent with the observations of RX J0822-4300 and IE 1207.4-5209. 
More rapid cooling for these two models of proton superfluidity is provided by the neutrino 
emission due to Cooper pairing of protons which occurs at t ~ 50-100 yr. 

Thick lines in the left panel of Fig. El demonstrate the additional effect of neutron pairing 
nsl in the crust. Comparing three thick lines, one can see that crustal neutron pairing 
noticeably accelerates only very slow cooling of low-mass neutron stars with strong proton 
pairing pi in their cores (Yakovlev et al. 2001, 2002). In that case the neutrino luminosity 
due to Cooper pairing of neutrons in the stellar crust at t < 3 x 10 5 yr may dominate the total 
neutrino luminosity of the stellar core. Moreover, at t > 3 x 10 5 yr crustal neutron pairing 
reduces the heat capacity of the crust. Both effects accelerate the cooling and decrease T s °°, 
violating the interpretation of the two hottest sources, RX J0822-4300 and IE 1207.4-5209. 
Any model of weaker crustal superfluidity will only bring cooling curves closer to thin ones 
and simplify the interpretation of the observations. 

On the other hand, for moderately strong (p2) or moderate (p3) proton pairing in the 
core, the effects of strong crustal neutron pairing on the cooling of middle-aged neutron stars 
(10 3 yr < t < 10 5 yr) are almost negligible. The neutrino emission due to crustal Cooper 
pairing of neutrons can noticeably accelerate the cooling and decrease T s °° only during the 
internal thermal relaxation stage (t < 100 yr). 

The right panel of Fig. El demonstrates that the observations of RX J0822-4300 and 
IE 1207.4-5209 can be explained by adopting any model of proton pairing (pi, p2 or p3), 
model nsl of crustal superfluidity, and the presence of an accreted envelope of the mass 
AM = 10~ 8 M (thin lines). 

Let us note, that the upper dot-and-short-dashed cooling curve goes higher than is needed 
to interpret the observations of the young and hottest source RX J0822-4300. Accordingly, 
following Yakovlev et al. (2002) (also see Potekhin et al. 2003), we may assume the presence 
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of a thinner accreted envelope (e.g., AM ~ 10 _11 M Q ) to interpret the observations of RX 
J0822-4300 and IE 1207.4-5209 (for the combination of pi and nsl pairing). The stronger 
proton core superfluidity, the less massive accreted envelope is needed for the interpretation 
of the data for these two stars. 

In order to explain the old and warmest sources, PSR B1055-52 and RX J0720. 4-3125, 
we will treat them as low-mass stars with the iron surface and proton pairing p2 in the core 
(or similar model of pairing with the peak of critical temperature T™ ax > 10 9 K). Moreover, 
the presence of any crustal neutron pairing (for example, nsl; thick solid lines in Fig. EJ), 
does not violate the interpretation of these sources. Note that proton pairing p3 (thick 
long-dashed lines) is less appropriate for the interpretation of these sources than pairing 
p2. Therefore, we adopt proton pairing p2 as the basic model for a new cooling scenario. 
Obviously, any model of stronger proton pairing (with higher T cp (p)) is better consistent 
with the observations. 

5 ACCRETED ENVELOPES AND COOLING OF NEUTRON STARS 

Figure El illustrates the effects of the accreted envelopes of the mass AM = 10~ 8 M on the 
cooling of neutron stars with different masses and the same nucleon pairing (models p2, 
ntl, and nsl). For comparison, we present also the cooling curves for stars with iron surface 
(thick solid lines) and the same nucleon superfluidity (also see the right panel of Fig. Q]). 
Note that the effect of crustal superfluidity on the cooling of such stars is unimportant. 

In the left panel of Fig. El we present our traditional cooling curves T s °°(t) and compare 
them with the data on the surface temperatures. On the right panel we show the temporal 
evolution of the surface thermal luminosity L^(t) and compare it with the data (Table 1). 
Both representations of the same cooling processes are seen to be in a reasonably good 
agreement although the data on are generally less certain and seem to be currently 
less conclusive (because, as a rule, the luminosity of the selected sources is determined less 
accurately than their surface temperature as discussed in Section EJ) . 

Figure El shows a strong rise of cooling curves for neutron stars with accreted envelopes 
at the neutrino cooling stage (t < 3 x 10 4 yr) and their steep decrease at the photon cooling 
stage. Their photon stage starts earlier than for stars with the iron surface. Assuming the 
presence of accreted envelopes, we can explain the observations of the young and hottest 
neutron stars, RX J0822-4300 and IE 1207.4-5209, treating them either as low-mass or 
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Figure 3. Cooling curves of neutron stars with different masses and nucleon pairing p2, ntl, and nsl (Fig. Q versus ob- 
servations. Thin dot-and-dashed (M = 1.111 Mq), long dashed (M = 1.37 Mq), and short dashed (M = 1.994 Mq) curves 
are calculated including the effects of accreted envelopes (AM = 10 — 8 M). Three thick solid curves are the same as the solid 
curves in the right panel of Fig.0 The left and right panels show the same cooling curves but as functions T s °°(t) and LJ°(t), 
respectively. 

as medium-mass stars. In contrast, the observations of the old and warmest objects, PSR 
B1055-52 and RX J0720. 4-3125, can be explained only by treating them as low-mass stars 
with the iron surfaces and with moderately strong (or strong) proton pairing inside. 

It was shown by Chang & Bildsten (2003, 2004) that the mass of light elements may 
decrease with time, particularly due to diffusive nuclear burning. The characteristic burning 
time r can be considered as an additional cooling regulator. Following Chang & Bildsten 
(2003, 2004) and Page et al. (2004) we assume that the mass of light elements decreases 
with time as AM(t) = AMq exp(— t/r), where AMq is the initial mass. 

Figure 0] illustrates the effect of variable mass of the accretion envelope on the cooling 
of the M = I.IIIMq neutron star. All cooling curves are calculated assuming nucleon 
pairing p2, ntl, and nsl. The thick solid line is our typical cooling curve for a low-mass 
superfluid star without any accreted envelope. We use two values of the initial mass of light 
elements AMq/M = 10 -9 and 10~ 7 and present thus three pairs of cooling curves for three 
characteristic times r. When r is lower than the time of the transition from the neutrino 
cooling stage to the photon stage (r < 3 x 10 4 yr) we obtain a smooth transition of the 
cooling track from the regime of highest temperatures in young stars to the regime of lower 
temperatures in old stars (cf. the curves for M = 1.111 M & in Figs. HHandEJ). This effect 
has been pointed out by Page et al. (2004). At r > 3 x 10 5 yr cooling curves merge into 
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Figure 4. Cooling curves of 1.111 Mq neutron stars with nucleon pairing p2, ntl, and nsl versus observations. The thick 
solid curve is the same as in Fig. El other curves are calculated assuming an exponential decay of accreted envelopes with 
characteristic times r = 0.4 kyr (short-dashed lines), 4 kyr (long-dashed lines), and 40 kyr (dot-and-dashed lines) for two 
values of the initial accreted mass, AMq/M = 10 — 9 and 10 — 7 . 



the limiting curve obtained for constant AM = AM . In the intermediate case of 3 x 10 4 yr 
< t < 3 x 10 5 yr the cooling curves gradually approach this limiting curve with the increase 
of r. As seen from Fig. EJ by assuming any r in the range 10 3 yr < r < 10 yr one can 
explain the observations of all neutron stars hottest for their age by one cooling curve. Note 
also that the value AMq/M = 10~ 9 is too small to explain the observations of young and 
hottest neutron stars, especially RX J0822-4300, at any r. 

As remarked by Chang & Bildsten (2004), an accreted envelope of a pulsar can become 
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thinner owing to the excavation of ions from the stellar surface by a pulsar wind at a rate 
M ~ 2f2 2 mj/i/ec, where Q is the pulsar spin frequency, \x is the magnetic moment, and m; 
is the ion mass. For an ordinary pulsar with the spin period ~ 0.1 s, fi ~ 10 30 G cm 3 , and 
helium surface we would have the surface mass loss AM ex ~ 6 x 1CT 12 M Q in t ~ 10 5 yr, too 
small to affect the cooling of a star with the initial helium layer of AM > 1O~ 1O M . For 
a pulsar with much higher magnetic field and/or faster rotation the effect may be stronger 
and affects the cooling. 



6 CONCLUSIONS 

We have extended the scenario of neutron star cooling proposed by Gusakov et al. (2004a) 
taking into account the effects of accreted envelopes and crustal singlet-state pairing of 
neutrons. As stressed in Section^ this scenario is different from the minimal cooling scenario 
of Page et al. (2004). 

The general idea of the minimal cooling scheme is that the enhanced neutrino emission, 
required for the interpretation of observation of neutron stars coldest for their age, is provided 
by the neutrino emission due to Cooper pairing of neutrons. In this case the direct Urea 
process or similar enhanced neutrino processes in kaon- condensed, pion-condensed, or quark 
matter can be forbidden in neutron stars of all masses. 

As in Gusakov et al. (2004a), the proposed cooling scenario imposes stringent constraints 
on the density dependence of the critical temperature T cn (p) for triplet-state neutron pairing 
in the stellar core. They result from the comparison of theoretical cooling curves with the 
data on the three most important "testing sources" , PSR J0205+6449, RX J0007.0+7303, 
and the Vela pulsar (Sect. EJ). By tuning our phenomenological model of triplet-state neu- 
tron pairing in the stellar core we obtain a noticeable dependence of the cooling on neutron 
star mass. It enables us to explain all the data by single combination of models for nucleon 
superfluidity. Assuming the presence of accreted envelopes we obtain two additional param- 
eters to regulate the cooling, which are the initial envelope mass AMo and its characteristic 
burning time r (Chang & Bildsten 2003). 

Our interpretation implies the presence of moderately strong proton pairing (T™ ax > 10 9 K) 
and moderate triplet-state neutron pairing (with T™f x ~ 6 x 10 8 K) in neutron star cores. 
Also, we have taken into account the effect of strong singlet-state neutron pairing (T™ ax ~ 

7 x 10 9 K) in the stellar crust. However, as shown in Sects. H] and the effect of crustal 
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superfluidity is unimportant for cooling middle-aged neutron stars with moderately strong 
proton pairing in their cores. 

We need proton superfluidity to explain the observations of the neutron stars hottest 
for their age. However, in contrast to the cooling scenario of Gusakov et al. (2004a), our 
new cooling scenario does not require too strong proton pairing. In fact, we can explain the 
observations of the old and warmest stars, PSR B1055-52 and RX J0720. 4-3125, by treating 
them as low-mass neutron stars (without accreted envelopes) with moderately strong proton 
pairing in their cores. Such phenomenological models for proton pairing are consistent with 
recent microscopic calculations of proton critical temperatures by Zuo et al. (2004) and 
Takatsuka & Tamagaki (2004). 

The young and hottest neutron stars, RX J0822-4300 and IE 1207.4-5209, can also be 
treated as low-mass stars with the same moderate proton superfluidity in their cores but 
assuming the presence of accreted envelopes. The smaller the mass of the envelope, required 
for the interpretation of these sources, the stronger proton pairing should be assumed. 

As discussed above, we need neutron pairing ntl (or similar) to explain the observations 
of the stars coldest for their age. However, as has been demonstrated by Gusakov et al. 
(2004b), cooling curves are not too sensitive to exchanging neutron and proton superfluidities 
(T cp (p) # T cn (p)) in neutron star cores. Therefore, we would also be able to explain the data 
in the scenario with moderately strong neutron and moderate proton pairing in the stellar 
cores. 

Neutron star cooling can also be affected by surface magnetic fields and by some reheating 
mechanisms in neutron star interiors. We have not discussed the effects of magnetic fields 
(although they are incorporated in our cooling code). The main reason is that these effects 
are weaker than the effects discussed above (for ordinary cooling isolated neutron stars of 
non-magnetar type; see, e.g., Yakovlev et al. 2002 for a detailed discussion of this point). 
Internal reheating mechanisms (see, e.g., Page 1998a,b, and references therein), for instance, 
the reheating due to the viscous dissipation of differential rotation, are relatively weak and 
model dependent; they become important at the photon cooling stage. No reheating is 
required to explain the data in our cooling scenario. What is more important, that most 
elaborated model equations of state of dense matter (Akmal, Pandharipande & Ravenhall 
1998) predict the operation of the direct Urea process in most massive stable neutron stars. 
This should lead to the existence of new classes of cooling neutron stars. The scenario with 
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the open direct Urea process (which can be called the extended minimal cooling scenario) 
has been studied by Gusakov et al. (2005). 

It is important that the same physics of neutron star interiors, which is tested by obser- 
vations of isolated (cooling) neutron stars, can also be tested by observations of accreting 
neutron stars in soft X-ray transients (e.g., Yakovlev, Levenfish & Haensel 2003) basing on 
the hypothesis of deep crustal heating of such stars (Brown, Bildsten & Rutledge 1998) by 
pycnonuclear reactions in accreted matter (Haensel & Zdunik 1990). The observations of 
soft X-ray transients in quiescent states indicate (Yakovlev, Levenfish & Gnedin 2005) the 
existence of rather cold neutron stars (first of all, SAX J1808. 4-3658) inconsistent with the 
model of neutron star structure proposed in the present paper. However, these observational 
indications are currently inconclusive (e.g., Yakovlev et al. 2005). If confirmed in future 
observations, they could give stronger evidence against the proposed scenario than new ob- 
servations of cooling neutron stars. In this case the extended minimal cooling scenario may 
appear to be more perspective. 
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